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We can understand many recently-discovered features of flavor evolution in dense, self-coupled 
supernova neutrino and antineutrino systems with a simple, physical scheme consisting of two 
quasi-static solutions. One solution closely resembles the conventional, adiabatic single neutrino 
Mikheyev-Smirnov-Wolfenstein (MSW) mechanism, in that neutrinos and antineutrinos remain in 
mass eigenstates as they evolve in flavor space. The other solution is analogous to the regular pre- 
cession of a gyroscopic pendulum in flavor space, and has been discussed extensively in recent works. 
Results of recent numerical studies are best explained with combinations of these solutions in the 
following general scenario: (1) Near the neutrino sphere, the MSW-like many-body solution obtains. 
(2) Depending on neutrino vacuum mixing parameters, luminosities, energy spectra, and the matter 
density profile, collective flavor transformation in the nutation mode develops and drives neutrinos 
away from the MSW-like evolution and toward regular precession. (3) Neutrino and antineutrino 
flavors roughly evolve according to the regular precession solution until neutrino densities are low. 
In the late stage of the precession solution, a stepwise swapping develops in the energy spectra of 
and fp/i/r. We also discuss some subtle points regarding adiabaticity in flavor transformation in 
dense neutrino systems. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

Although it was well known that neutrino-neutrino for- 
ward scattering could generate sizable flavor refractive in- 
dices in dense- neutrino environments [E 0, B S S @| i it 
was thought generally that the large matter density near 
the supernova neutrino sphere would suppress neutrino 
self-coupling effects. However, it recently became appar- 
ent that neutrino self-coupling could alter drastically the 
way neutrinos and antineutrinos evolve in flavor space in 
the supernova environment 0, H, , even for the small 
neutrino mass-squared differences inferred from experi- 
ments 

Assuming coherent neutrino flavor transformation and 
the efficacy of a mean field approach [l^, [llj, and us- 
ing the magnetic spin analogy for 2x2 neutrino oscilla- 
tions and the technique of the "corotating frame" (dis- 
cussed later in the text), Ref. [l^l was able to show that 
neutrinos can experience collective flavor transformation 
even in the presence of a dominant matter field. The 
first large-scale numerical simulations of neutrino fiavor 
transformation in the late-time supernova environment 
with correlated neutrino trajectories were discussed in 
Refs. 13, 14]. The results of these simulations clearly 
revealed the collective nature of neutrino oscillations in 
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supernovae. These effects are not explainable by the con- 
ventional Mikheyev-Smirnov-Wolfenstein (MSW) effect 

The numerical simulations [H, [l^ showed that neu- 
trinos of all energies and propagating along all trajec- 
tories can experience rapid flavor oscillations on short 
time/distance scales in the late-time supernova hot- 
bubble region. This type of collective flavor transfor- 
mation was actually first discussed in the context of the 
early universe [H, [l^, . Ref. [IH showed that a uni- 
form mono-energetic neutrino gas initially in pure Ve and 
De states is equivalent to a gyroscopic pendulum in fla- 
vor space, and that the rapid collective neutrino flavor 
transformation in supernovae and in the early universe 
corresponds to the nutation of this flavor pendulum. 

The simulations discussed in Refs. [13, [3] also revealed 
a stunning feature of collective neutrino flavor transfor- 
mation: i^e's appear to swap their energy spectra with 
i^fj,/r 's at energies below or above (depending on the neu- 
trino mass hierarchy) a transition energy Eq ■ It was sug- 
gested in Ref. [l3l | that this stepwise swapping in neutrino 
energy spectra is related to a collective mode in which the 
relative phases of the two flavor components of neutrinos 
oscillate at the same rate. Ref. [22!| argued that this col- 
lective mode is related to the regular precession of a fla- 
vor pendulum. This was done by solving the equations 
for regular precession mode of a gas of mono-energetic 
neutrinos initially in pure j/g and Dp states. Following 
the suggestion proposed in Ref. Ref. [13] further 

elucidates why a stepwise swapping in neutrino energy 
spectra can result from neutrinos evolving in a collec- 
tive precession mode. Ref. found a way to obtain 
the regular precession solution for neutrino gases with 
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any give n neutrino number density and energy spectra. 
Ref. [23] also pointed out that the transition energy Eq 
was determined by the conserved lepton number in the 
regular precession solution. This solution was termed an 
"adiabatic solution" in Ref. ,23]. However, we will show 
later in this paper that the flavor evolution in the regular 
precession solution is actually not adiabatic. 

Inspired by the method that Ref. [23] used to find the 
regular precession solution, we show that a quasi-static 
solution for given electron density and neutrino number 
density can be derived in a similar way. This solution 
is very much like the conventional adiabatic MSW solu- 
tion, except that the background neutrinos contribute an 
additional refractive index to a given propagating neu- 
trino. We calculate this MSW-like solution as well as 
the regular precession solution under the same condi- 
tions as the single-angle numerical simulations discussed 
in Refs. [l^, [l^. The comparison between these solu- 
tions and the numerical results renders a simple picture 
for neutrino flavor transformation in supernovae: Neutri- 
nos initially follow the MSW-like solution near the neu- 
trino sphere before being driven away from this solution 
by collective flavor transformation in the nutation mode. 
Subsequently, neutrinos roughly follow the regular pre- 
cession solution with some nutation. As the neutrino 
and antineutrino densities decrease, stepwise swapping 
appears in the neutrino energy spectra. 

This paper is organized as follows. In Sec. |lT] we de- 
rive the MSW-like solution and discuss the conditions 
under which flavor evolution in this scenario is adiabatic. 
In Sec. mil we briefly recapitulate the derivation of the 
regular precession solution and discuss the effects of the 
matter field on this solution. In Sec. IIVI we compare the 
MSW-like solution and the regular precession solution to 
the numerical results and discuss the origin of stepwise 
swapping in neutrino energy spectra. We also elaborate 
on the adiabaticity of collective neutrino flavor evolution. 
In Sec.|V]we give our conclusions. 



The corresponding NFIS in the flavor basis is 
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The flavor basis wavefunction for an antineutrino is de- 
flned as 



(3) 



where ap^ and a^^ are the amplitudes for the antineutrino 
to be De and Dt, respectively. The corresponding NFIS 
in the flavor basis is 



-, / 2Re(ap^a* ) 
2Im(a,,a;j 
2 \la,j2-la,j2 



(4) 



As one will see later, the special definition of the flavor 
basis wavefunction and NFIS for the antineutrino allows 
one to write the e.o.m. for both neutrinos and antineu- 
trinos in a unified manner. 

The fiavor evolution of a neutrino is represented by the 
precession of the corresponding NFIS around an effective 
external field H. The e.o.m. of NFIS can be written 
in a form similar to that for a magnetic spin: 



dt 



Si X H,: 



(5) 



where subscript i denotes the initial physical state (fla- 
vor species and momentum) , or "mode" , of the neu- 
trino/ antineutrino. As Si precesses around H;, the pro- 
jection of Si on oscillates. This represents variation in 
the flavor content of the neutrino. Here 
are the unit vectors in the flavor basis. 

In the absence of ordinary matter and other neutrinos 
(in vacuum), the effective fleld experienced by NFIS Sj is 



and 



II. MSW-LIKE SOLUTION FOR DENSE 
NEUTRINO GASES 



Hi — ujill\ 



(6) 



where the vacuum fleld is 



A. Equations of motion 



-e!, sin26'v -t- e{ cos26'v, 



(7) 



Here we briefly recapitulate a discussion of the equa- 
tions of motion (e.o.m.) for neutrino flavor transfor- 
mation. We use the notation of neutrino flavor isospin 
(NFIS) (see Ref. [l2| for detailed discussions). A NFIS is 
defined as the expectation value of the Pauli spin opera- 
tor of the flavor wavefunction for a neutrino or antineu- 
trino. The flavor basis wavefunction for a neutrino in the 
2x2 mixing scheme is 



(1) 



where a^^ and a^^ are the amplitudes for the neutrino 
to be in v,, and another flavor state, say v-r, respectively. 



and LOi is 



OJi 



±- 



2^ 



(8) 



In this notation 6^ G (0, tt/2) is the effective 2x2 vacuum 
mixing angle, Sm^ > is the mass-squared difference 
between two neutrino vacuum mass eigenstates, and Ei 
is the energy of the neutrino. The plus sign in Eq. ([5]) 
is for neutrinos, and the minus sign is for antineutrinos. 
Throughout this paper we will assume that sin26'v = 
sin 26*13 ^1- In this convention one has 0^ ^ 1 for 
the normal neutrino mass hierarchy and 9^ ~ tt/2 for 
the inverted one. Eqs. ([5]) and ^ for NFIS can be 
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pictured as the precession of around Hy with angular 
velocity uji. 

In the presence of ordinary matter, the forward scat- 
tering of neutrinos on electrons induces different refrac- 
tive indices for neutrinos of electron flavor from those 
for other flavors. As a result, in the presence of matter, 
NFIS Si will also tend to precess around the vector field 

He EE -V2GFneel. (9) 

Here Gf is the Fermi constant and Ue is the net electron 
number density. 

For dense neutrino gases, forward neutrino-neutrino 
scattering couples the flavor evolution of different neu- 
trino modes. With all the contributions taken into ac- 
count, the effective field experienced by NFIS is 



Hi = CJiHv + He -f ^ TljSj, 



where 



(10) 



(11) 



is the coupling coefhcient of two neutrino modes for 
isotropic neutrino gases, and Uj is the population (num- 
ber density) of neutrino mode j. We define the normal- 
ized distribution function as 

/.-^. (12) 



where 



(13) 



Using Eq. ((TO)) we can express Eq. ([14]) in the explicit 
component form in the vacuum mass basis: 

Si.x = -;^[He.x + f^iySx), (16a) 
Si.z ^ TTTri^i + ^e.z + fJ'i^Sz). (16c) 

Here Sx(y^z) are the components of the total NFIS 

S = ^n,s, = n*°'^/,s, (17) 



in the same basis, and 



Hi = {He.x + IJ-uSxY + ifJ'iySy)'^ + {uJi + He,z + Mi/'S'z)^- 

(18) 

Summing Eq. (|16|) over index i with weight Ui = n^^^fi 
we obtain 



— —^{He,x + fJ-uSx) ~77~' 



tot 

tot 
v 

tot 



(19a) 
(19b) 



Sz 



nT eJ, 



Y^^{u;, + He,z+f^uSz). (19c) 



is the total neutrino number density. 



Eqs. (fT9a|) and (figb]) imply Sy = if He ^ 0. In this 
case, and for any given n*°* and Ue, one can solve for Sx 
and Sz using the equations 



B. MSW-like solution 

In the NFIS notation, the instantaneous light (heavy) 
mass eigenstate of a neutrino is represented by a NFIS 
completely aligned (anti-aligned) with its effective field. 
(Note that the effective energy of NFIS s,; is £°*^ = 
—Si • H,;.) If the flavor evolution of neutrinos is fully 
adiabatic, a neutrino initially in the light (heavy) mass 
eigenstate will stay in the same instantaneous mass eigen- 
state, and a NFIS initially aligned (anti-aligned) with its 
effective field will stay aligned (anti-aligned) with the ef- 
fective field. In such a limit one has 



H^ 



(14) 



where the alignment factor = +1 (—1) if NFIS Si is 

aligned (anti-aligned) with its effective field H^. 

It is convenient to work in the vacuum mass basis in 

f 



which the unit vectors ^x(y,z) are related to e 

= cos 20V + sin 26'v, 
f 



Hv 



■Ay,-) by 

(15a) 
(15b) 

-e^ sin 26lv + cos 26'v. (15c) 



Sx — 2 {-He,x ~l" fJ'i'Sx) 



E 



V iHe,x + f^i^Sx)'^ + {uJi + He,z + t^-i^SzY 



(20a) 



fv-, I''-"' 



2 ^ yj {He,x + iJ-uSxY + (^i + He.z + Mi/S'^)^ 

(20b) 

The components of NFIS s^ can then be obtained from 
Eq. p6)) . We note that the counterpart of Eq. ((20|) in the 
flavor basis was discussed in connection with adiabatic 
flavor transformation of supernova neutrinos in Ref . [2^ . 

We also note that the vanishing of the z-component of 
S in the flavor basis, i.e.. 



si = Sz cos 26lv -I- Sx sin 26lv = 0, 



(21) 



corresponds to an MSW-like resonance at which the to- 
tal NFIS represents a maximally mixed state. Using 
Eqs. (|20ap and (j20bp we rewrite the above equation as 

(E 2^?v + f E #1 - 0' (22) 
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where H: 



-^/2G-pne- In the Hmit n*; 



oo, 



all 



NFIS's are either aligned or anti-aligned with the total 
NFIS: 



2S 



(note that /i^ < 0). In the same limit, Eq. 
gives a normalization condition 



(23) 

20a[) simply 



2 ~ ^ 



e,)h - 1, (24) 



and Eq. ([22|) reduces to 



Wsync cos 26'v + Hi ., ~ 0, 



(25) 



where 



'^sync 



2 

^ w.in^Si • S ~ ^(-e,)/,tj, (26) 



^tot 

25* 



where i^i is the angle between and Hy in the MSW- 
like solution. Eq. ([^ is the condition for the MSW-like 
solution to be adiabatic. 

We note that the MSW-like solution becomes more adi- 
abatic with larger neutrino densities. This is because Hi 
is the energy gap between the instantaneous light and 
heavy mass eigenstates of neutrino mode i, and it in- 
creases with n|;°*. 

We also note that Eq. (|28)) only gives the necessary con- 
dition for the neutrino system to adiabatically follow the 
MSW-like solution. This is because in deriving Eq. ([28|) 
we have assumed that is described by the MSW-like 
solution in the first place. This is true for the conven- 



tional MSW solution where ; 



= 0. If 7 



is large, how- 



ever, all NFIS's can move in a collective manner, even in 
the presence of the matter field |12i|. As a result, the to- 
tal NFIS S and the effective field Hi = tJ^Hv -I- He -f ^^S 
may not follow the MSW-like solution at all. We will 
further elaborate on this point in Sec. IIVDI 



is the synchronization frequency [255. Eq. (|25|) implies 
that all neutrinos go through an MSW-like resonance 
at approximately the radius where a neutrino with en- 
ergy E = 5w? /2\ujsyac\ would for a conventional, single- 
neutrino MSW resonance. 

Eqs. (|19ap and (jl9bp are not independent ii He = 0. 
Assuming S is static with given n|f * , one can still choose 
a new coordinate system in which Sy vanishes. Naively, 
it might seem practical to solve for Sx and Sz using 
Eqs. (pna|l and (|20bp . However, from Eqs. ^ and HI 
one can show that the lepton number 



H. 



2Sz 



(27) 



is conserved for Hg = [2l| . This is true even if is 
a function of time [23|. Therefore, and are over 
constrained by Eqs. (|20ap . (|20b|) and ((77)1 . In this case, 
a new unknown is expected. In fact, a regular precession 
solution exists in the absence of the matter field. In such 
a solution, S precesses around Hy and the precession 
angular velocity ujpr is the new unknown variable to be 
found (see Sec. IIIIIl. 



Adiabatic condition 



The MSW-like solution inherent in Eqs. (pOal) and 
(|20b[) is a quasi-static solution for any given Tig and n*°' . 
In this solution, as Hi changes with varying Ue and n*°* , 
Si can lag behind. In other words, as the ensemble of 
neutrinos evolves, misalignment between Si and Hi can 
develop. Consequently, Si will tend to precess around Hi 
with angular speed Hi. Therefore, in order for the fla- 
vor evolution of neutrino mode i to follow the MSW-like 
solution, one must have 



= Hr^ 



dt 



« 1, 



(28) 



III. REGULAR PRECESSION SOLUTION FOR 
DENSE NEUTRINO GASES 

A. Regular precession solution 

In the absence of a matter field (iJg — 0), the e.o.m. of 
NFIS's possess cylindrical symmetry around Hy. A so- 
lution with the same symmetry is expected to exist for 
any given n*°*. In this solution, which we term the "regu- 
lar precession solution", all NFIS's must precess steadily 
around Hy and 



— Sj = Sj X WprHy. 



(29) 



For a gas of mono-energetic neutrinos initially in pure 
and Ve states, this represents the regular precession of 
the gyroscopic pendulum in flavor space [22| . 
Combining Eqs. HU]) and we have 

S, X [{Ui - Wpr)Hy + /vS] = 0, (30) 

which means that Si is either aligned or anti-aligned with 



H, 



(31) 



We note that Hi is the effective field experienced by NFIS 
Si in a frame whose i 
frame according to 



Si in a frame whose unit vectors e^^^^ ^-j rotate in the static 



d 2v iv 



(32) 



Therefore, the regular precession solution can be derived 
by means similar to those used to get the MSW-like so- 
lution ^23]. 

Eq. ^ shows that the effective fields for aU NFIS's 
and, therefore, all NFIS's themselves, reside in a static 
plane in the corotating frame which is spanned by Hy 
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and S. Without loss of generality, we assume this coro- 
tating plane to be the e^-e^ plane, and write out the 
components of NFIS s^: 



Ci UJi — Wpr + l-liySz 



(33a) 
(33b) 



where subscripts x and z indicate the projections of the 
vectors onto e^. and e^, respectively, and the alignment 
factor is = +1 (—1) for aligned (anti-aligned) with 
Hi. We sum Eq. ()33a|) with weight n^°^fi and obtain 



1 = 



E 



(34) 



Summing Eq. (|33bp with weight 'n}°*' fi and using 
Eq. ([341), we obtain 



E 



(Wi - UJp^ + + {\luSxf 



(35) 

For any given values of total neutrino number density 
n*°* and lepton number £, Eq. (P7)) specifies 5*2. One 
can then obtain Sx and o^p^ from Eqs. (|34p and (|35p . 
The components of each NFIS in the corotating frame 
are determined by Eq. 



B. Pseudo-regular precession condition 

Because the angle between and Hy varies with n*°* , 
it is not possible for to stay in the regular precession 
solution. In this case, will tend to precess around 
H; in the corotating frame. This precession corresponds 
to nutation in the static frame. Ref. (2^ has argued 
that the regular precession solution can be an excellent 
approximation to the actual evolution of the system if 
the inverse of the nutation time scale is much larger than 
|di9i/di|, where di is the angle between and Hy in 
the regular precession solution. As is parallel to H^, 
•Qi is equivalent to the angle between Hj and Hy for 
the solution. If this is the case, the amplitude of the 
nutation will be small. This behavior is similar to the 
pseudo-regular precession of a gyroscope. 

We note that the time scale of the nutation of in the 
static frame is the same as the period of the precession 
of Si around H^ in the corotating frame, which is Hi. 
Therefore, the condition for NFIS to be in the pseudo- 
regular precession is 



1% 



< 1. 



(36) 



We note that Eq. (136]) only gives the necessary condi- 
tion for neutrino systems to follow the regular precession 



solution. This is because in deriving Eq. (|36p we have 
assumed that H^ follows the regular precession solution 
and rotates much slower than does . This can be true if 
the nutation of all NFIS's are not correlated. In collective 
flavor transformation, however, the motion of the total 
NFIS S and, therefore. Hi is correlated with individual 
NFIS Si. As a result. Hi moves at a rate comparable 
to that of Si, and the flavor evolution may not exactly 
follow the regular precession solution. 

The pseudo-regular precession condition in Eq. (j36|) 
was first proposed in Ref. [23] as an "adiabatic condi- 
tion" . The notion "adiabatic condition" here can be con- 
fusing or even misleading. Because the effective field for 
a NFIS in the static frame can be different from that 
in a corotating frame, the alignment/anti-alignment of a 
NFIS with its effective field in one frame does not guar- 
antee the alignment/anti-alignment in the other frame. 
We note that the flavor evolution of a neutrino is con- 
sidered to be adiabatic if, e.g., it stays in the light mass 
eigenstate or, equivalently, the corresponding NFIS stays 
aligned with its effective field in the static frame. In con- 
trast, Eq. is a necessary condition for all NFIS's to 
remain aligned or anti-aligned with their effective fields 
in the appropriate corotating frame. As we will show at 
the end of Sec. IIVDI the flavor evolution of neutrinos is 
not adiabatic in the conventional sense if they stay in the 
regular precession mode, which satisfies the condition in 
Eq. 



C. Effect of the matter field 

Ref. HI] argues that, so long as sin26'v is small, ordi- 
nary matter has no effect other than changing the preces- 
sion angular velocity Wpr of the system. This is because, 
if sini^i 3> sin26'v, NFIS Si essentially sees He as parallel 
to Hy, and 



H, 



CJiHy 

o;^Hy 



where 



(37a) 
(37b) 



(38) 



Eq. (j37bp takes a form similar to that in the e.o.m. of 
NFIS Si in the absence of the matter field. We note that 



where 



"pr 



(39) 



and Wpi- is the precession velocity in the absence of the 
matter field. Therefore, Eq. (|30p would hold and the 
precession solution would obtain in the corotating frame 
with angular velocity if H^^x were indeed ignorable. 

In the case of regular precession, Si would be aligned or 
anti-aligned with its effective field Hi in the frame which 
rotates around Hy with angular velocity lo'^,. In this 
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corotating frame, the NFIS would experience an addi- 
tional field He which has magnitude |-ffe,2:| and precesses 



around Hy with angular velocity 



As long as 



« K.\ = 



(40) 



He would only introduce small perturbation in the mo- 
tion of Si and the regular precession solution is a good ap- 
proximation. The regular precession approximation fails 
if 



-'pr 



0, 



(41) 



which approximately describes the conventional MSW 
resonance condition for a neutrino with energy E — 
6m? /2\u)p^\. 

In addition to the above effect of the matter field on 
the collective precession in the regime of high neutrino 
number density, the matter field can also cause individual 
neutrino modes to deviate from the collective precession 
mode at their MSW resonances if the neutrino number 
density is low. This can be understood as follows. The 
effective field H^ for NFIS in the static frame can be 
written as the sum of two fields: 



H 



MSW 



(42) 



and /ii/S^, where Sj_ is the component of the total NFIS 
S perpendicular to e^. As n|;°* and rie decrease, Hmsw 
will rotate, and, consequently, neutrino mode i will en- 
counter an MSW-like resonance when 



0. 



(43) 



At the same time, precesses around Hy with angu- 
lar velocity Wpr. If the rotation speed of Hmsw is slow 
enough for the MSW resonance to be adiabatic while at 
the same time | /Ziyis'^ < Wpr, NFIS will follow Hmsw 
through the resonance and, therefore, will deviate from 
the collective precession mode. Obviously, the lepton 
number C is not conserved in this process. 



A. Supernova model 

In the single-angle simulations discussed in Refs. [H, 
IT3 |. the flavor evolution of all neutrinos is assumed to 
be the same as that of the neutrinos propagating along 
radial trajectories. Under this assumption, the neutrino- 
neutrino coupling strength is taken to be = —2^/2G-p 
as in isotropic neutrino gases. The neutrino-neutrino in- 
tersection angle dependence in the current-current weak 
interaction is partially taken into account by introducing 
an effective number density of neutrinos [1, 5]. The ef- 
fective number density rfj^ of neutrinos of energy E at 
radius r is defined as 



_ D{r/R,)L,UE) 



(44) 



In Eq. ((44|) , Ri, is the radius of the neutrino sphere (here 
taken to be = 11 km), the geometric factor 



D{rlR.) ^ - 




(45) 



incorporates the geometric coupling and dilution of 
anisotropic neutrino beams as a function of radius r, 
and {E^) are the luminosity and average energy of the 
neutrino species, respectively, and fu{E) is the (normal- 
ized) energy distribution function for neutrinos. 

In the single-angle simulations, luminosities for all neu- 
trino species are taken to be the same, and two typi- 
cal late-time supernova neutrino luminosity values L^, = 
5 X 10^^ erg/s and 10^^ erg/s have been used. The energy 
distribution function fu{E) for neutrinos is taken to be 
of the Fermi-Dirac form with two parameters (T,^,riv), 



ME) 



1 



1 



E^ 



F2(7;,)T3exp(£;/T,-r;,) + r 



(46) 



where we take the degeneracy parameter to be t;,^ = 3, 
is the neutrino temperature, and 



IV. NEUTRINO FLAVOR TRANSFORMATION 
IN SUPERNOVAE 

The problem of neutrino flavor transformation under 
realistic supernova conditions is very difficult to solve. 
This is largely a result of the anisotropic nature of the 
neutrino and antineutrino distribution functions in the 
supernova environment. However, it has been demon- 
strated numerically that most of the qualitative features 
of the supernova neutrino oscillation problem are cap- 
tured by the so-called "single-angle" simulations (see be- 
low) [ij, [3, H^. In this section we will show explic- 
itly that the results from single-angle simulations can be 
explained by the combination of MSW-like solution and 
the regular precession solution discussed in the preceding 
sections. 



X Ax 



exp(a; — rj) -\- \ 



(47) 



The values of Ti, for various neutrino species are deter- 
mined from our chosen average energies: {E^^) = 11 
MeV, {E^:) = 16 MeV, and {E,^) = (S^J = \e,^) = 
{E^^) = 25 MeV Q. 

The simulations use a simple density profile. The net 
electron density at radii sufficiently above the neutrino 
sphere is taken to be @ 



Mns ruN 



'PI 



1-4 -3 



(48) 



where Yg is the electron fraction (here taken to be Ye — 
0.4), cjs — 11/2 is the statistical weight in relativistic par- 
ticles, A^NS = 1-4M0 is the mass of the neutron star, rnjv 
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is the mass of the nucleon, mpi is the Plank mass, and 
S = 140 is the entropy in units of Boltzmann's constant 
per baryon. 

The matter density near the neutrino sphere is much 
larger than what Eg. (1481) calculates. The simulations 
discussed in Refs. [ij, |lJ] have adopted an exponential 
density profile in the region near the neutrino sphere. We 
note that a large rig will keep neutrinos in their initial 
flavor eigenstates. For the MSW-like solution, the 
profile in Eq. is large enough to keep neutrinos in 
their initial flavor states to r ^ R^, , and the exponential 
density profile near the neutrino sphere will not affect 
our analysis below. 



Comparison of numerical and analytical 
solutions 



For 2x2 flavor mixing, neutrinos starting in pure Ve 
(Pe) and Vr (i^t) states with the same energy evolve in 
the same way and can be viewed as the same neutrino 
mode. This is because NFIS and — follow the same 
e.o.m. As all neutrinos are assumed to be in pure flavor 
eigenstates at the neutrino sphere, a neutrino mode i is 
uniquely designated by uji [defined in Eq. ([5])], and we 
can take Vl- f dto. 

We define the total neutrino number density at radius 
r as 



nT (r)^ I < {E,r)~ < {E,r)\ dE 
Jo 

+ / \nf^{E,r)-nf^{E,r)\dE. 
Jo 



Note that nj;°* defined above takes advantage of the 
equivalent flavor evolution of (f'e) and u-r (ut), and 
is not the simple sum of number densities of all neutri- 
nos. We define the NFIS mode distribution function as 



~ ^\dE/M^ 



\nl«{E,R,)^nl^{E,R,)\ ii lu > 0, 
\nf{E,R,)-nf{E,R,)\ ii uj < 0. 

(50) 

Because /i^ < and i^g's are dominant in number in 
supernovae, the vector ii^S is initially in the direction 
opposite to e^. Taking this into account, we have 



He + fii,S 

H,,, 



(51) 



at the neutrino sphere, where is the total effective 
field for neutrino mode uj. Eq. (|5ip is true for any neu- 
trino mode UJ. 

Noting that ly'e/i^T and De/i'T are represented by NFIS's 
aligned and anti-aligned with , respectively, we find the 
alignment factor for neutrino mode uj to be 

^ r -sgn(nf (E, R,) - nf{E, R,)) if > 0, 
1 sgninf(E, R,) - nfiE, R,)) if w < 0, 



where sgn(^) = ^/|^| is the sign of ^. 

Using Eqs. (gH), ^ and ^ we are able to find 
the MSW-like solution in the same supernova model as 
adopted in the single-angle simulations. Taking 5m? — 
3 X 10"'^ eV^ (close to the value associated with atmo- 
spheric neutrino oscillations), 9^ = 0.1 (normal neutrino 
mass hierarchy) or 7r/2 — 0.1 (inverted neutrino mass hi- 
erarchy), = 5 X 10^^ erg/s, we compute the average 
NFIS 



n 



S_ 

tot 



(53) 



in the MSW-like solution. We plot {s±) = \{sx) \ and (s^), 
the perpendicular and z components, respectively, of the 
average NFIS in the vacuum mass basis, as dot-dashed 
lines in Fig. [1] 

We also obtain the regular precession solution under 
the same conditions. For the inverted mass hierarchy 
case, we determine the value of (s^) in the regular pre- 
cession solution using the lepton number C = 2{sz) = 
2Sz/n^°^ at the neutrino sphere [see discussions around 
Eq. (|54p ]. For the normal mass hierarchy case, we deter- 
mine the value of (sz) using the transition energy Eq [see 
discussions around Eq. ([6T|) ]. Note that the swapping of 
energy spectra above/below Eq (stepwise swapping) is a 
natural result of the regular precession solution. With 
^2 = nl°^Sz we solve Eqs. dH]) and for Sx (or S± in 
the static frame) and ujp^. The values of {s±) = S±/n}°^ 
and (sz) corresponding to this solution for both the in- 
verted and normal mass hierarchies are plotted as solid 
(49) lines in Fig. [T] 

For comparison with the MSW-like and regular preces- 
sion solutions, we extract the values of {s±) and (s^) in 
the single-angle simulations presented in Ref. [l^ . These 
numerical results are plotted as dashed lines in Fig. [TJ 

For the inverted mass hierarchy case [Fig. [Tfa-b)] we 
observe that {s±) in the numerical simulation follows the 
MSW-like solution until radius rx — 88 km, but then 
abruptly jumps out of the MSW-like solution. There- 
after, {s±) roughly follows the regular precession solu- 
tion. The reason for this change in the behavior of the 
NFIS's is discussed extensively in Ref. [22] . Briefly speak- 
ing, supernova neutrinos naturally form a "bipolar sys- 
tem" with the neutrinos forming two groups of NFIS's 
pointing in opposite directions in flavor space [31 . This 
bipolar system is roughly equivalent to a gyroscopic pen- 
dulum in flavor space |21i] . For the inverted mass hier- 
archy case, the gyroscopic pendulum is initially in the 
highest position, and (n*°*)~^ plays the role of (variable) 
gravity. For a given value of internal spin (determined 
by J'e-j^e asymmetry), the gyroscopic pendulum will not 
fall until the "gravity" is large enough and ri*°* drops be- 
low some critical value n^. If the gravity increases very 
slowly, so that the symmetry around Hy is preserved, 
the gyroscopic pendulum will precess steadily. This cor- 
responds to the regular precession solution. However, 
the symmetry is usually broken near the point where 
(52) n*°*(r) ~ n^. Below this critical density, the pendu- 
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FIG. 1: (Color online) Plots of {s±{r)) (left panels) and {sz{r)} (right panels) for both the inverted (upper panels) and normal 
(lower panels) neutrino mass hierarchies. The dashed lines show numerical simulation results from Ref. The solid lines 

show the regular precession solution. The dot-dashed lines show the MSW-like solution. The luminosity for all neutrino species 
is taken to be L„ = 5 x lO^'^ erg/s. 



lum will develop significant nutation in addition to the 
precession. 

We note that the fast oscillations of {s±) around the 
regular precession solution correspond to the nutation 
mode. We also note that (sz) roughly stays constant in 
the region r < 115 km and is 



(54b) 

~ -0.1. (54c) 

This is because both the collective precession and nuta- 
tion modes conserve the lepton number C = 2(sz). The 
value of {sz) is subsequently changed because of matter 
effects, which do not conserve C. However, (sj_) in the 
numerical solution still roughly follows that of the regu- 
lar precession solution, even after (sz) starts to evolve. 
This seems to imply that the neutrino system is still de- 



scribed by the regular precession solution, except for the 
neutrino modes which have dropped out of the collective 
precession. 

For the normal mass hierarchy case [Fig. [BJc-d)] we 
observe that (sj_) and (sz) in the numerical simulation 
follow the MSW-like solution until radius r ~ 83 km. At 
this point, (s) has almost completely flipped its direction. 
This corresponds to the complete flavor transformation 
in both Vf. ^ Vt and C'e ^ Dt channels. In the flavor pen- 
dulum analogy, the gyroscopic pendulum has been raised 
from the lowest position to the highest position, and one 
expects it to fall when n*°*(r) < nl- Indeed, for large 
enough radii, both {sj_) and (sz) leave the MSW-like so- 
lution and start to oscillate around the regular precession 
solution with (sz) — —0.07. 

We have also obtained the three kinds of solutions 
for a smaller neutrino luminosity, L^, = 10^^ erg/s, with 
other parameters unchanged. These results are plotted 
in Fig. O Comparing Fig. [T]and[21 one can see that the 
small luminosity cases are very similar to their large lu- 
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FIG. 2: (Color online) Same as Fig. [T] except that the luminosity for all neutrino species is taken to be Li, = lO^'^ erg/s in this 
figure. 



minosity counterparts. However, we note that, for the 
inverted mass hierarchy case, (sj_) leaves the MSW-hke 
solution at rx — 63 km, which is smaller than the rx 
value in the large luminosity case. This is because with 
smaller L^, the total neutrino number density reaches 
earlier. For the normal mass hierarchy case, our cal- 
culation procedure cannot find the MSW-like solution in 
a small radius interval immediately beyond r ~ 91 km. 
We have taken (s^) ~ —0.16 for the regular precession 
solution in this case. We also note that the value of 
(sz) changes significantly in both small luminosity cases. 
This is because with smaller neutrino number densities, 
the matter field drives more neutrinos or antineutrinos 
off the track of the collective precession solution. 



C. Stepwise swapping in neutrino energy spectra 



aligned with effective fields in the appropriate corotating 
frame (see Sec. IIII A]) . It was first suggested in Ref. [IjI 
and further explained in Ref. [23 | that, if this align- 
ment in the corotating frame is maintained, Ve and I'r 
will swap their energy spectra for energy smaller (larger) 
than a transition energy Eq in the normal (inverted) 
mass hierarchy case. This phenomenon is also referred 
to as a "spectral split" and the transition energy is de- 
termined by the conservation of lepton number £ [23j . 
This stepwise swapping of neutrino energy spectra has 
been observed in the numerical simulations presented in 
Refs. [ll,[ll. 

The phenomenon of the stepwise spectrum swapping 
(or spectral split) can be understood as follows. In the 
regular precession solution, NFIS s^^ stays aligned or anti- 
aligned with and s„ — e^'H.^/2H^. The effective field 
in the corotating frame becomes 



As we have seen in Sec. IIVBI supernova neutrinos 
evolve initially according to the MSW-like solution, but 
subsequently roughly follow the regular precession solu- 
tion. In the latter solution, all NFIS's are aligned or anti- 



= (c^ - )H^ 



(55) 



as v}"^ 0, where uj^^. is the value of ujpi- at = 0. 
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Therefore, one has 



(56) 



In other words, the spectrum of e^jScj.z is spht into two 
parts when njf* — > 0: the lower part {uj < uj^^) takes 
the value —1/2 while the upper part {uj > ujp^) takes the 
value +1/2. 

For sin 20V ^ 1 the NFIS for mode oj in the limit of 
n*°* = is 

(si^)„t_ot^o = Hv(sw,z)„tot^Q, (57a) 
- sgn(w - w° Je^^ , (57b) 

where S = +1 (—1) for the normal (inverted) mass hi- 
erarchy. Eq. (j57bp means that supernova neutrinos are 
almost in pure flavor eigenstates when n|^°* — > 0. We 
note that neutrinos are in flavor eigenstates at the neu- 
trino sphere and 



(Si^)r=fl:„ 



(58) 



Therefore, the probability for a neutrino mode lu to stay 
in its original state is 



P^i,{uj) = - + 2(s^)r=_R„ • (Sc^)„tot^o, 

- I [l-2sgn(cj-LjJ^,)] . 



(59a) 
(59b) 



According to Eq. (j59bp . a stepwise swapping can occur 
in the final neutrino energy spectra. For lu^j. > and 
9v — 7r/2 (inverted mass hierarchy), Ve and Vr will swap 
their energy spectra for energies above 



En — 



2<- 



(60) 



For tjpj. > and ^ 1 (normal mass hierarchy) , V(. and 
Vr will swap their energy spectra for energies below Eq. 

Assuming stepwise spectrum swapping, the value of 
Wpj. can be determined from (sz) without finding the com- 
plete regular precession solution ^2^: 



2(.s.) = 



(61) 



Using the values of (s^) in the regular precession solutions 
[Eq. ([54])]. we obtain Ec — 8.4 MeV for the inverted mass 
hierarchy cases. This value agrees well with the numer- 
ical simulation results (see Fig. 9(c) in Ref. 13]). For 
the normal mass hierarchy scenarios, we actually obtain 
(sz) in the regular precession solution by using Eq. (pT|) 
and the values of Eq in the numerical results (Fig. 9(a) 
in Ref. [l3|) which are approximately 7.8 MeV and 9.5 
MeV for the large and small neutrino luminosity cases, 
respectively. As shown in Figs. andlD^c), the pre- 
cession solutions determined by these values of Eq agree 



well with the numerical simulation results. For the pa- 
rameters we have chosen, the stepwise swapping occurs 
only in the neutrino sector. This is because Ve is the 
dominant (in terms of number flux) neutrino species in 
supernovae. 



D. Adiabaticity 

The flavor evolution of neutrino mode lo is adiabatic so 
long as the angle between and H^^ remains constant. It 
can be shown that the flavor evolution of neutrino mode 
UJ is adiabatic if H^^ rotates at a rate much slower than 
H^, or, in other words. 



H,„ X 



dH„ 



dt 



(62) 



(see, e.g., Ref. fl^l). Adiabaticity is guaranteed if the 
condition in Eq. (|62[) is satisfied. However, when neutrino 
densities are large, the system can evolve adiabatically 
even if the condition in Eq. ([62]) is not satisfied. 

If neutrinos experience collective flavor transformation, 
the total NFIS S, and therefore, effective fleld Huj ~ 
LuHy + He + /iyS move in correlation with the motion 
of individual NFIS In this case, can move very 
fast, with staying aligned or anti- aligned with IH^, 
and the flavor evolution is, therefore, still adiabatic. We 
expect this to be the case as neutrinos transition from 
the MSW-like solution through the collective nutation 
mode to the regular precession solution. As a result, the 
alignment factor for neutrino mode uj remains constant 
during the transition. This was implicitly assumed in the 
discussion on stepwise spectrum swapping in Sec. IIV CI 

On the other hand, we note that flavor evolution 
through the regular precession solution is not adia- 
batic. Although the pseudo-regular precession condition 
in Eq. (|36p is in a form similar to the adiabatic condition 
in Eq. ((62)) . it simply guarantees that NFIS s^^ remains 
aligned or anti-aligned with the effective field in the 
corotating frame, but not with in the static frame. 
Therefore, neutrino mode uj does not stay in the same 
mass eigenstate throughout its evolution. One may con- 
sider two neutrino modes uj^^ + Auj and uj^^. — Auj which 
are initially in pure states. If the flavor evolution is 
adiabatic for both neutrinos, both NFIS's should stay 
anti-aligned with their effective fields in the static frame 
throughout the evolution. However, we know that this is 
not the case. Because of the stepwise swapping in neu- 
trino energy spectra, one of the two NFIS's will become 
aligned with its effective field in the static frame when 
n*°* 0. Therefore, the flavor evolution of neutrinos in 
the regular precession solution cannot be adiabatic for all 
neutrinos. 
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V. CONCLUSIONS 

We have derived a pair of equations from which it is 
possible to find a quasi-static MSW-Uke solution for an 
isotropic neutrino gas with specified electron and neu- 
trino number densities. This solution is the natural ex- 
tension of the conventional MSW solution, but includes 
neutrino self-coupling. We have shown that the condi- 
tion for this MSW-like solution to be adiabatic is a nec- 
essary condition for neutrinos to follow such a solution. 
We have also discussed the previously discovered regular 
precession solution and the scenarios where it may break 
down in the presence of ordinary matter. 

We have compared the results of the detailed simula- 
tions in Refs. [13, with the corresponding MSW-like 
and regular precession solutions. This comparison clearly 
shows that many features in the numerical simulations, 
including the stepwise swapping in neutrino energy spec- 
tra, can be explained by combinations of these two ana- 
lytical solutions. We have also discussed the adiabaticity 
of flavor evolution in dense neutrino gases. 

We emphasize that the discovery of the MSW-like so- 
lution and the regular precession solution by no means 
obviates the need for further numerical simulations. For 
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